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Abstract 

Heavy Ion Collisions (HIC) represent a unique tool to probe the in-medium nuclear interaction 
in regions away from saturation. In this report we present a selection of new reaction observables 
in dissipative collisions particularly sensitive to the symmetry term of the nuclear Equation of 
State (I so — EoS). We will first discuss the Isospin Equilibration Dynamics. At low energies 
this manifests via the recently observed Dynamical Dipole Radiation, due to a collective neutron- 
proton oscillation with the symmetry term acting as a restoring force. At higher beam energies 
Iso-EoS effects will be seen in Imbalance Ratio Measurements, in particular from the correlations 
with the total kinetic energy loss. For fragmentation reactions in central events we suggest to 
look at the coupling between isospin distillation and radial flow. In Neck Fragmentation reactions 
important I so — EoS information can be obtained from the correlation between isospin content 
and alignement. The high density symmetry term can be probed from isospin effects on heavy 
ion reactions at relativistic energies (few AGeV range). Rather isospin sensitive observables are 
proposed from nucleon/cluster emissions, collective flows and meson production. The possibility 
to shed light on the controversial neutron/proton effective mass splitting in asymmetric matter is 
also suggested. A large symmetry repulsion at high baryon density will also lead to an "earlier" 
hadron-deconfinement transition in n-rich matter. A suitable treatment of the isovector interaction 
in the partonic EoS appears very relevant. 
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1 Introduction: The Elusive Symmetry Term of the EoS 



The symmetry energy E sym appears in the energy density e(p, p 3 ) = e(p) + pE sym (p 3 / p) 2 + 0(p 3 /p) 4 + 
.., expressed in terms of total (p = p p + p n ) and isospin (p 3 = p p — p n ) densities. The symmetry 
term gets a kinetic contribution directly from basic Pauli correlations and a potential part from the 
highly controversial isospin dependence of the effective interactions. Both at sub-saturation and supra- 
saturation densities, predictions based of the existing many-body techniques diverge rather widely, see 



We recall that the knowledge of the EoS of asymmetric matter is very important at low densities 
( e.g. neutron skins, pigmy resonances, nuclear structure at the drip lines, neutron distillation in 
fragmentation, neutron star formation and crust) as well as at high densities ( e.g. neutron star mass- 
radius relation, cooling, hybrid structure, transition to a deconfined phase, formation of black holes). 
Several observables which are sensitive to the Iso-EoS and testable experimentally, have been suggested 
[SI SI El El [8] . We take advantage of new opportunities in theory (development of rather reliable 
microscopic transport codes for HIC) and in experiments (availability of very asymmetric radioactive 
beams, improved possibility of measuring event-by-event correlations) to present new results that are 
constraining the existing effective interaction models. We will discuss dissipative collisions in a wide 
range of beam energies, from just above the Coulomb barrier up to the AGeV range. Isospin effects 
on the chiral/deconfmement transition at high baryon density will be also discussed. Low to Fermi 
energies will bring information on the symmetry term around (below) normal density, while intermediate 
energies will probe high density regions. The transport codes are based on mean field theories, with 
correlations included via hard nucleon-nucleon elastic and inelastic collisions and via stochastic forces, 
selfconsistently evaluated from the mean phase-space trajectory, see (5j [9j [US [EE]- Stochasticity is 
essential in order to get distributions as well as to allow the growth of dynamical instabilities. 

Relativistic collisions are described via a fully covariant transport approach, related to an effective 
field exchange model, where the relevant degrees of freedom of the nuclear dynamics are accounted for 
[SI E21 EH EEU E3, EE] • We will have a propagation of particles suitably dressed by self-energies that will 
influence collective flows and in medium nucleon-nucleon inelastic cross sections. The construction of 
an Hadron — EoS at high baryon and isospin densities will finally allow the possibility of developing a 
model of a hadron-deconfinement transition at high density for an asymmetric matter [17] . The problem 
of a correct treatment of the isospin in a effective partonic EOS will be stressed. 

We will always test the sensitivity of our simulation results to different choices of the density and 
momentum dependence of the Isovector part of the Equation of State {I so— EoS). In the non-relativistic 
frame the potential part of the symmetry energy, C(p), [5]: 



is tested by employing two different density parametrizations, Isovector Equation of State (Iso-Eos) 
[31 [TH], of the mean field: i) ^ = 482 - 1638p, (MeVfm 3 ), for "Asysoft" EoS: E sym /A(pot) has 



a weak density dependence close to the saturation, with an almost flat behavior below po and even 
decreasing at suprasaturation; ii) a constant coefficient, C = 32MeV, for the "Asystiff" EoS choice: 
the interaction part of the symmetry term displays a linear dependence with the density, i.e. with a 
faster decrease at lower densities and much stiffer above saturation. The isoscalar section of the EoS is 
the same in both cases, fixed requiring that the saturation properties of symmetric nuclear matter with 
a compressibility around 220MeV are reproduced. 



[l,2j. 
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2 Isospin Equilibration 



2.1 The Prompt Dipole 7-Ray Emission 

The possibility of an entrance channel bremsstrahlung dipole radiation due to an initial different N/Z 
distribution was suggested at the beginning of the nineties [HI [20]. After several experimental evi- 
dences, in fusion as well as in deep-inelastic reactions, [211 [22] and refs. therein, we have now a good 
understanding of the process and stimulating new perspectives from the use of radioactive beams. 

During the charge equilibration process taking place in the first stages of dissipative reactions be- 
tween colliding ions with different N/Z ratios, a large amplitude dipole collective motion develops in 
the composite dinuclear system, the so-called Dynamical Dipole mode. This collective dipole gives rise 
to a prompt 7-ray emission which depends: i) on the absolute value of the intial dipole moment 

NZ 

D{t = 0) = — - \R z (t = 0) - R N (t = 0)| = 
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being Rz = and Rn = S '^"^ the center of mass of protons and of neutrons respectively, while 

Rp and Rp are the projectile and target radii; ii) on the fusion/deep-inelastic dynamics; iii) on the 
symmetry term, below saturation, that is acting as a restoring force. 

A detailed description is obtained in mean field transport approaches, [23, 24J. We can follow the time 



evolution of the dipole moment in the r-space, D(t) = ^-(Rz—Rn) and in p— space, DK(t) = — 
with P p (P n ) center of mass in momentum space for protons (neutrons), just the canonically conjugate 
momentum of the D(t) coordinate, i.e. as operators [D(t), DK(t)] = %h. A nice "spiral-correlation" 
clearly denotes the collective nature of the mode, see Fig.l. 

We can directly apply a bremsstrahlung approach, to the dipole evolution given from the Landau- 
Vlasov transport [24] , to estimate the "prompt" photon emission probability {E 1 = Two): 

dP 2e 2 

^ = 37^ |/y ' H|2 ' (3) 

where D"(uj) is the Fourier transform of the dipole acceleration D"(t). We remark that in this way it 
is possible to evaluate, in absolute values, the corresponding pre-equilibrium photon emission. 

We must add a couple of comments of interest for the experimental selection of the Dynamical 
Dipole: i) The centroid is always shifted to lower energies (large deformation of the dinucleus); ii) A 
clear angular anisotropy should be present since the prompt mode has a definite axis of oscillation 
(on the reaction plane) at variance with the statistical GDR. In a recent experiment the prompt 
dipole radiation has been investigated with a 4ir gamma detector. A strong dipole-like photon angular 
distribution (6> 7 ) = W [l + a 2 P2{cos9~ / )], # 7 being the angle between the emitted photon and the beam 
axis, has been observed, with the a 2 parameter close to —1, see [22J. 

At higher beam energies we expect a decrease of the direct dipole radiation for two main reasons 
both due to the increasing importance of hard NN collisions: i) a larger fast nucleon emission that will 
equilibrate the isospin during the dipole oscillation; ii) a larger damping of the collective mode due to 
np collisions. 

The use of unstable neutron rich projectiles would largely increase the effect, due to the possibility 
of larger entrance channel asymmetries [22]. In order to suggest proposals for the new RIB facility 
Spiral 2, [26] we have studied fusion events in the reaction 132 Sn + 58 Ni at 10 AMeV, [251 12Z]- We 
espect a Monster Dynamical Dipole, the initial dipole moment D(t = 0) being of the order of 50fm, 
about two times the largest values probed so far, allowing a detailed study of the symmetry potential, 
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Figure 1: Dipole Dynamics at lOAMeV, b = Afm centrality. Left Panels: Exotic "132" system. Upper: 
Time evolution of dipole moment D(t) in real space; Lower: Dipole phase-space correlation (see text). 
Right Panels: same as before for the stable "124" system. Solid lines correspond to Asysoft EoS, the 
dashed to Asystiff EoS. 



I32 Sn + 58 Ni lOAMeV b=4fm l24 Sn + 58 Ni lOAMeV b=4fm 




Figure 2: Left Panel, Exotic "132" system. Power spectra of the dipole acceleration at b = 4fm (in 
c 2 units). Right Panel: Corresponding results for the stable "124" system. Solid lines correspond to 
Asysoft EoS, the dashed to Asystiff EoS. 



below saturation, responsible of the restoring force of the dipole oscillation and even of the damping, 
via the fast neutron emission. 

In Figure 1 we report some global informations concerning the dipole mode in entrance channel. 
In the Left-Upper panel we have the time evolution of the dipole moment D(t) for the "132" system 
at b = Afm. We notice the large amplitude of the first oscillation but also the delayed dynamics for 
the Asystiff EOS related to a weaker isovector restoring force. The phase space correlation (spiraling) 
between D(t) and DK(t), is reported in Fig.l (Left-Lower). It nicely points out a collective behavior 
which initiates very early, with a dipole moment still close to the touching configuration value reported 
above. This can be explained by the fast formation of a well developed neck mean field which sustains 
the collective dipole oscillation in the dinuclear configuration. 

The role of a large charge asymmetry between the two colliding nuclei can be seen from Fig.l (Right 
Panels), where we show the analogous dipole phase space trajectories for the stable 124 SYz+ 58 Ni system 
at the same value of impact parameter and energy. A clear reduction of the collective behavior is 
evidenced. 

In Fig. 2 (Left Panel) we report the power spectrum, | D"(u) | 2 in semicentral "132" reactions, 
for different I so — EoS choices. The gamma multiplicity is simply related to it, see Eq.Q. The 
corresponding results for the stable "124" system are drawn in the Right Panel. As expected from 
the larger initial charge asymmetry, we clearly see an increase of the Prompt Dipole Emission for the 



4 



exotic n-rich beam. Such entrance channel effect will be enhanced, allowing a better observation of the 
Iso-EoS dependence. In fact from Fig. 2 we see also other isospin effects. 

A detailed analysis can be performed just using a simple damped oscillator model for the dipole 
moment D(t) = D(to)e t( - UJO+l ^ T,>t , where D(to) is the value at the onset of the collective dinuclear response, 
lvq the frequency and r the damping rate. The power spectrum of the dipole acceleration is given by 

I D>) |»= M+j/iW (4) 

(u - uj ) 2 + 1/r 

which from Eq.flS]) leads to a total yield proportional to u t(uJq + l/T 2 )D(t ) 2 ~ uJqtD(1; ) 2 since 
ujqt > 1. We remind that in the Asystiff case we have a weaker restoring force for the dynamical dipole 
in the dilute "neck" region, where the symmetry energy is smaller [5] . This is reflected in lower values of 
the centroids as well as in reduced total yields, as shown in Fig. 2. The sensitivity of ujq to the stiffness 
of the symmetry energy will be amplified by the increase of D(to) when we use exotic, more asymmetric 
beams. 

The prompt dipole radiation angular distribution is the result of the interplay between the collective 
oscillation life-time and the dinuclear rotation. In this sense we expect also a sensitivity to the I so — Eos 
of the anusotropy, in particular for high spin event selections, [28] . 

In the Asysoft choice we expect also larger widths of the "resonance" due to the larger fast neutron 
emission. We note the opposite effect of the Asy-stiffness on neutron vs proton emissions. The latter 
point is important even for the possibility of an independent test just measuring the N/Z of the pre- 
equilibrium nucleon emission, [29J . 

2.2 Isospin Equilibration at the Fermi Energies 

In this energy range the doorway state mechanism of the Dynamical Dipole will disappear and so we 
can study a direct isospin transport in binary events. This can be discussed in a compact way by means 
of the chemical potentials for protons and neutrons as a function of density p and isospin / [30]. The 
p/n currents can be expressed as 

h/n = D p p/r y P -D I p/n VI (5) 

with Dp, the drift, and D p , n the diffusion coefficients for transport, which are given explicitely in ref. 
[30] . Of interest here are the differences of currents between protons and neutrons which have a simple 
relation to the density dependence of the symmetry energy 

D p -D p oc 4I dEsym 

Dl-Dl ex ApE sym . (6) 

Thus the isospin transport due to density gradients, i.e. isospin migration, depends on the slope of the 
symmetry energy, or the symmetry pressure, while the transport due to isospin concentration gradients, 
i.e. isospin diffusion, depends on the absolute value of the symmetry energy. 

We can discuss the asymmetries of the various parts of the reaction system (gas, PLF/TLF's , and 
in the case of ternary events, IMF's ). In particular, we study the so-called Imbalance Ratio [31], which 
is defined as 

2(x M -x^) 

>T ~ p _ X L) > W 

with x eq = \{x H + x L ). Here, x is an isospin sensitive quantity that has to be investigated with respect 
to equilibration. In this work we consider primarily the asymmetry {3 = (N — Z) / (N + Z) , but also 
other quantities, such as isoscaling coefficients, ratios of production of light fragments, etc, can be of 
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Figure 3: Left Panel. Imbalance ratios for Sn + Sn collisions for incident energies of 50 (left) and 35 
AMeV (right) as a function of the impact parameter. Signatures of the curves: iso-EoS stiff (solid 
lines), soft (dashed lines); MD interaction (circles), MI interaction (squares); projectile rapidity ( full 
symbols, upper curves ), target rapidity ( open symbols, lower curves ). Right Panel. Imbalance ratios 
as a function of relative energy loss for both beam energies. Upper: Separately for stiff (solid) and soft 
(dashed) iso-EoS, and for MD (circles and squares) and MI (diamonds and triangles) interactions, in 
the projectile region (full symbols) and the target region (open symbols). Lower: Quadratic fit to all 
points for the stiff (solid), resp. soft (dashed) iso-EoS. 



interest [7]. The indices H and L refer to the symmetric reaction between the heavy (n-rich) and the 
light (ra-poor) systems, while M refers to the mixed reaction. P, T denote the rapidity region, in which 
this quantity is measured, in particular the PLF and TLF rapidity regions. Clearly, this ratio is ±1 in 
the projectile and target regions, respectively, for complete transparency, and oppositely for complete 
rebound, while it is zero for complete equilibration. 

In a simple model we can show that the imbalance ratio mainly depends on two quantities: the 
strength of the symmetry energy and the interaction time between the two reaction partners. Let 
us take, for instance, the asymmetry (3 of the PLF (or TLF) as the quantity x. As a first order 
approximation, in the mixed reaction this quantity relaxes towards its complete equilibration value, 
(3 eg = ((3h + Pl)/2 as 

(3™ T = p e i + (p H ' L - (3 £ i) e ~ t/T , (8) 

where t is the time elapsed while the reaction partners are interacting (interaction time) and the damping 
t is mainly connected to the strength of the symmetry energy. Inserting this expression into Eq. (f7|) , 
one obtains Rp T = ±e - */ r for the PLF and TLF regions, respectively. Hence the imbalance ratio can 
be considered as a good observable to trace back the strength of the symmetry energy from the reaction 
dynamics provided a suitable selection of the interaction time is performed. 

The centrality dependence of the Imbalance Ratio, for (Sn,Sn) collisions, has been investigated in 
experiments [32] as well as in theory [30, 33J. We propose here a new analysis which appears experi- 
mentally more selective. The interaction time certainly influences the amount of isospin equilibration, 
see Eq. (6) and refs. [3U1 EH] • Longer interaction times should be correlated to a larger dissipation. It 
is then natural to look at the correlation between the imbalance ratio and the total kinetic energy loss. 
In this way we can also better disentangle dynamical effects of the isoscalar and isovector part of the 
EoS, see [M]. 

It is seen in Fig. 3 that the curves for the asy-soft EoS (dashed) are generally lower in the projectile 
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Figure 4: Left Panel. The N/Z of the liquid (left) and of the gas (right) phase is displayed as a function 
of the system initial N/Z. Full lines and symbols refer to the asy stiff parameterization. Dashed lines 
and open symbols are for asysoft. Right Panel. The fragment N/Z (see text) as a function of the kinetic 
energy. Left: Asy stiff; Right: Asysoft. 



region (and oppositely for the target region), i.e. show more equilibration, than those for the asy-stiff 
EoS. In order to emphasize this trend we have, in the lower panel of the figure, collected together all 
the values for the stiff (circles) and the soft (squares) iso-EoS, and fitted them by a quadratic curve. It 
is seen that this fit gives a good representation of the trend of the results. 

The difference between the curves for the stiff and soft iso-EOS in the lower panel then isolates the 
influence of the iso-EoS from kinematical effects associated with the interaction time. It is seen, that 
there is a systematic effect of the symmetry energy of the order of about 20 percent, which should be 
measurable. The correlation suggested in Fig.3 should represent a general feature of isospin diffusion, 
and it would be of great interest to verify experimentally. 



3 Isospin Distillation with Radial Flow 

In central collisions at 30-50 MeV/A, where the full disassembly of the system into many fragments is 
observed, one can study specifically properties of liquid-gas phase transitions occurring in asymmetric 
matter [35J, [361 E3 [HI E] . For instance, in neutron- rich matter, phase co-existence leads to a different 
asymmetry in the liquid and gaseous phase: fragments (liquid) appear more symmetric with respect 
to the initial matter, while light particles (gas) are more neutron-rich. The amplitude of this effect 
depends on specific properties of the isovector part of the nuclear interaction, namely on the value and 
the derivative of the symmetry energy at low density. 

This investigation is interesting in a more general context: In heavy ion collisions the dilute phase 
appears during the expansion of the interacting matter. Thus we study effects of the coupling of 
expansion, fragmentation and distillation in a two-component (neutron-proton) system [38J . 

We focus on central collisions, 6 = 2 fin, considering symmetric reactions between systems hav- 
ing three different initial asymmetry: n2 Sn + 112 Sn, 124 Sn + 124 Sn, 132 Sn + 132 Sn, with (N/Z) in = 
1.24,1.48,1.64, respectively. The considered beam energy is 50 MeV/A. 1200 events have been run for 
each reaction and for each of the two symmetry energies adopted (asy-soft and asystiff, see before) 
[38] . The average N/Z of emitted nucleons (gas phase) and Intermediate Mass Fragments (IMF) is 
presented in Fig. 4 (Left Panel) as a function of the initial (N/Z) in of the three colliding systems. One 
observes a clear Isospin-Distillation effect, i.e. the gas phase (right) more neutron-rich than the IMF's 
(left). This is particular evident in the Asysoft case due to the larger value of the symmetry energy at 
low density [5]. 

In central collisions, after the initial collisional shock, the system expands and breaks up into many 
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Figure 5: Correlation between N/Z of IMF and alignement in ternary events of the l24 Sn+ 64 Ni reac- 
tion at 35 AMeV. Left Panel. Exp. results: points correspond to fast formed IMFs (Viola-violation 
selection); histogram for all IMFs at mid-rapidity (including statistical emissions). Right Panel. 
Simulation results: squares, asysoft; circles, asy stiff. 

pieces, due to the development of volume (spinodal) and surface instabilities. The formation of a 
bubble- like configuration is observed, where the initial fragments are located |39j . 

Fragmentation originates from the break-up of a composite source that expands with a given velocity 
field. Since neutrons and protons experience different forces, one may expect a different radial flow for 
the two species. In this case, the N/Z composition of the source would not be uniform, but would 
depend on the radial distance from the center or mass or, equivalently, on the local velocity. This trend 
should then be reflected in a clear correlation between isospin content and kinetic energy of the formed 
IMF's, [38]. 

This observable is plotted in Fig.4 (Right Panel) for the three reactions. The behaviour observed 
is rather sensitive to the Iso-EoS. For the proton-rich system, the N/Z decreases with the fragment 
kinetic energy, expecially in the Asystiff case (left), where the symmetry energy is relatively small at 
low density. In this case, the Coulomb repulsion pushes the protons towards the surface of the system. 
Hence, more symmetric fragments acquire larger velocity. The decreasing trend is less pronounced in the 
Asysoft case (right) because Coulomb effects on protons are counterbalanced by the larger attraction 
of the symmetry potential. In systems with higher initial asymmetry, the decreasing trend is inversed, 
due to the larger neutron repulsion in neutron-rich systems. 

In conclusion, this analysis reveals the existence of significant, EoS-dependent correlations between 
the N/Z and the kinetic energy of IMF's produced in central collisions. 

4 Isospin Dynamics in Neck Fragmentation at Fermi Energies 

It is now quite well established that the largest part of the reaction cross section for dissipative collisions 
at Fermi energies goes through the Neck Fragmentation channel, with IMFs directly produced in the 
interacting zone in semiperipheral collisions on very short time scales [101 E]- We can predict interesting 
isospin transport effects for this new fragmentation mechanism since clusters are formed still in a 
dilute asymmetric matter but always in contact with the regions of the projectile-like and target-like 
remnants almost at normal densities. As discussed in Sect. 2. 2 in presence of density gradients the 
isospin transport is mainly ruled by drift coefficients and so we expect a larger neutron flow to the 
neck clusters for a stiffer symmetry energy around saturation, 0, E]. The isospin dynamics can be 
directly extracted from correlations between N/Z, alignement and emission times of the IMFs. The 
alignment between PLF — IMF and PLF — TLF directions represents a very convincing evidence 
of the dynamical origin of the mid-rapidity fragments produced on short time scales |42j. The form 



8 




Figure 6: 132Sn+124Sn at 400AMeV, central coll.Isospin content of nucleon and light ion emissions vs 
p t (upper) and kinetic energy (lower). Upper Panel: Asysoft; Lower Panel: Asy stiff. 



of the $ p zane distributions (centroid and width) can give a direct information on the fragmentation 
mechanism [43J. Recent calculations confirm that the light fragments are emitted first, a general feature 
expected for that rupture mechanism [33]. The same conclusion can be derived from direct emission 
time measurements based on deviations from Viola systematics observed in event-by-event velocity 
correlations between IMFs and the PLF/TLF residues [321 1321 133]- We can figure out a continuous 
transition from fast produced fragments via neck instabilities to clusters formed in a dynamical fission 
of the projectile(target) residues up to the evaporated ones (statistical fission). Along this line it would 
be even possible to disentangle the effects of volume and shape instabilities. A neutron enrichment of 
the overlap ("neck") region is expected, due to the neutron migration from higher (spectator) to lower 
(neck) density regions, directly related to the slope of the symmetry energy [45J. 

A very nice new analysis has been performed on the Sn + Ni data at 35 AMeV by the Chimera 
Collab.44.6J, see FigJ3J left panel. A strong correlation between neutron enrichement and alignement 
(when the short emission time selection is enforced) is seen, that can be reproduced only with a stiff 
behavior of the symmetry energy, Fig. 3 right panel (for primary fragments) [47] . This represents a clear 
evidence in favor of a relatively large slope (symmetry pressure) around saturation. We note a recent 
confirmation from structure data, i.e. from monopole resonances in Sn-isotopes |48j. 



5 Isospin Effects at High Baryon Density: Effective Mass 
Splitting and Collective Flows 

The problem of Momentum Dependence in the Isovector channel (I so — MD) is still very controversial 
and it would be extremely important to get more definite experimental information, see the recent refs. 
[491 13U1 13T] . Exotic Beams at intermediate energies are of interest in order to have high momentum 
particles and to test regions of high baryon (isoscalar) and isospin (isovector) density during the reaction 
dynamics. Our transport code has been implemented with a BGBD — like [341 [33] mean field with 
a different (n,p) momentum dependence, see [50] I3T] . This will allow to follow the dynamical effect 
of opposite n/p effective mass splitting while keeping the same density dependence of the symmetry 
energy [54] . 

We present here some preliminary results for reactions induced by 132 Sn beams on 12A Sn targets 
at 400^471^6^ [52]. For central collisions in the interacting zone we can reach baryon densities about 
1.7 — 1.8po m a transient time of the order of 15-20 fm/c. The system is quickly expanding and the 
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Figure 7: Transverse momentum dependence of the difference between proton and neutron V 2 flows, 
at mid-rapidity, in a semi-central reaction Au+Au at 400AMeV. 

Freeze-Out time is around 50fm/c. 

In Fig. 6 we present the (n/p) and 3 H/ 3 He yield ratios at freeze-out, for two choices of Asy-stiffness 
and mass splitting, vs. transverse momentum (upper curves) and kinetic energy (lower curves). In 
this way we can separate particle emissions from sources at different densities. We note two interesting 
features: i) the curves are crossing at p t ~ projectile — 2.13/m -1 ; ii) the effect is not much dependent 
on the stiffness of the symmetry term. The crossing nicely corresponds to a source at baryon density 
p ~ 1.6poj [50l[5T]. These data seem to be suitable to disentangle Iso — MD effects. 

Collective flows are very good candidates since they are expected to be very sensitive to the mo- 
mentum dependence of the mean field, see The transverse flow, Vi(y,p t ) = (^), provides 
information on the anisotropy of nucleon emission on the reaction plane. Very important for the reac- 

2 2 

tion dynamics is the elliptic flow, V2(y,pt) = ( Px p z Pv ). The sign of V 2 indicates the azimuthal anisotropy 
of emission: on the reaction plane (V 2 > 0) or out-of-plane (squeeze — out, V 2 < 0) [53]. We have 
then tested the Iso — MD of the fields just evaluating the Dif ference of neutron/proton transverse 
and elliptic flows Vi^,~ p \y , p t ) = V{ l 2 (y,p t ) — V p 2 (y,p t ) at various rapidities and transverse momenta 
in semicentral (b/b max = 0.5) 197 Au + 197 Au collisons at 400 AMeV . For the nucleon elliptic flows the 
mass splitting effect is evident at all rapidities, and nicely increasing at larger rapidities and transverse 
momenta, with more neutron flow when m* < m * From Fig. 7 we clearly see how at mid-rapidity the 
mass splitting effects are more evident for higher tranverse momentum selections, i.e. for high density 
sources. In particular the elliptic flow difference becomes negative when m* < m*, revealing a faster 
neutron emission and so more neutron squeeze out (more spectator shadowing). In correspondance the 
proton flow is more negative (more proton squeeeze out) when m* > m* It is however difficult to draw 
definite conclusions only from proton data. The measurement of n/p flow differences appears essential. 
Due to the difficulties in measuring neutrons, our suggestion is to measure the difference between light 
isobar flows, like 3 if vs. 3 He and so on. We expect to still see effective mass splitting effects. 



6 Isospin Effects on Meson Production in Relativistic Heavy 
Ion Collisions 

The phenomenology of isospin effects on heavy ion reactions at intermediate energies (few AGeV range) 
is extremely rich and can allow a "direct" study of the covariant structure of the isovector interaction 
in a high density hadron medium. We work within a relativistic transport frame, beyond the cascade 
picture, consistently derived from effective Lagrangians, where isospin effects are accounted for in the 
mean field and collision terms. We show that rather sensitive observables are provided by the pion/kaon 
production K°/K + yields). Relevant non-equilibrium effects are stressed. 

An effective Lagrangian approach to the hadron interacting system is extended to the isospin degree 
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of freedom: within the same frame equilibrium properties (EoS, [56]) and transport dynamics can be 
consistently derived. Within a covariant picture of the nuclear mean field, for the description of the 
symmetry energy at saturation (a 4 parameter of the Weizsaecker mass formula) (a) only the Lorentz 
vector p mesonic field, and (b) both, the vector p (repulsive) and scalar 5 (attractive) effective fields are 
included. In the latter case a rather intuitive form of the Symmetry Energy can be obtained [T2l 



1j F 



sym . 2 



f P fs [ Ef 



mr x 2 



PB- ■ (9) 



The competition between scalar and vector fields leads to a stiffer symmetry term at high density [T2jl5]. 

We present here observable effects in the dynamics of heavy ion collisions. We focus our attention 
on the isospin content of meson production. The starting point is a simple phenomenological version 
of the Non-Linear (with respect to the iso-scalar, Lorentz scalar a field) effective nucleon-boson field 
theory, the Quantum-Hadro-Dynamics [22] • According to this picture the presence of the hadronic 
medium leads to effective masses and momenta M* = M + S s , k*^ = fc M — S M , with E s , scalar 
and vector self-energies. For asymmetric matter the self-energies are different for protons and neutrons, 
depending on the isovector meson contributions. We will call the corresponding models as NLp and 
NLp5, respectively, and just NL the case without isovector interactions. For the more general NLpS 
case the self-energies of protons and neutrons read: 

£ s (p,n) = -fa<7(Ps) ± fsPs3, 

^(p,n) = f u) f T f P J^ (10) 

(upper signs for neutrons), where p s = p sp + p S n, j a = jp+jn>Ps3 = Psp-p sn , 3z = 3p ~ 3n are tlie total 
and isospin scalar densities and currents and f a ,w,p,s are the coupling constants of the various mesonic 
fields. cr(p s ) is the solution of the non linear equation for the a field [T21 15]. From the form of the scalar 
self-energies we note that in n-rich environment the neutron effective masses are definitely below the 
proton ones. 

For the description of heavy ion collisions we solve the covariant transport equation of the Boltz- 
mann type within the Relativistic Landau Vlasov (RLV) method, using phase-space Gaussian test 
particles [21] , and applying a Monte-Carlo procedure for the hard hadron collisions. The collision term 
includes elastic and inelastic processes involving the production/absorption of the A(1232MeU) and 
N*(1440MeV) resonances as well as their decays into pion channels, [T5] . 

Kaon production has been proven to be a reliable observable for the high density EoS in the isoscalar 
sector [581 EH] - Here we show that the K 0,+ production (in particular the K°/K + yield ratio) can be 
also used to probe the isovector part of the EoS, |TB1 ISO] . 

Using our RMF transport approach we analyze pion and kaon production in central 197 Au + 197 Au 
collisions in the 0.8 — 1.8 AGeV beam energy range, comparing models giving the same "soft" EoS for 
symmetric matter and with different effective field choices for E sym . 

Fig. [8] reports the temporal evolution of A ±,0 ' ++ resonances, pions (7r =l= ' ) and kaons (K +, °) for 
central Au+Au collisions at lAGeV. 

It is clear that, while the pion yield freezes out at times of the order of 50/m/c, i.e. at the final 
stage of the reaction (and at low densities), kaon production occurs within the very early (compression) 
stage, and the yield saturates at around 20 fm/c. From Fig. [8] we see that the pion results are weakly 
dependent on the isospin part of the nuclear mean field. However, a slight increase (decrease) in the 
7r~ (7r + ) multiplicity is observed when going from the NL to the NLp and then to the NLpS model, 
i.e. increasing the vector contribution f p in the isovector channel. This trend is more pronounced for 
kaons, see the right panel, due to the high density selection of the source and the proximity to the 
production threshold. Consistently, as shown in the insert, larger effects are expected for early emitted 
kaons, reflecting the early N/Z of the system. 
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Figure 8: Time evolution of the A ±,0,++ resonances and pions 7r ±,C) (left), and kaons (K +, ° (right) for a 
central (b = fm impact parameter) Au+Au collision at 1 AGeV incident energy. Transport calculation 
using the NL, NLp, NLpS and DDF models for the iso- vector part of the nuclear EoS are shown. The 
inset contains the differential K°/K + ratio as a function of the kaon emission time. 

When isovector fields are included the symmetry potential energy in neutron-rich matter is repulsive 
for neutrons and attractive for protons. In a HIC this leads to a fast, pre-equilibrium, emission of 
neutrons. Such a mean field mechanism, often referred to as isospin fractionation [5], is responsible 
for a reduction of the neutron to proton ratio during the high density phase, with direct consequences 
on particle production in inelastic NN collisions. Threshold effects represent a more subtle point. The 
energy conservation in a hadron collision in general has to be formulated in terms of the canonical 
momenta, i.e. for a reaction 1 + 2 — > 3 + 4 as s in = (fcf + k^) 2 = (k^ + k^) 2 = s out . Since hadrons are 
propagating with effective (kinetic) momenta and masses, an equivalent relation should be formulated 
starting from the effective in- medium quantities k*^ = k^ 1 — S M and m* = m + E s , where S s and S M are 
the scalar and vector self-energies, Eqs. (fT0|) . The self-energy contributions will influence the particle 
production at the level of thresholds as well as of the phase space available in the final channel. In 
fact the threshold effect is dominant and consequently the results are nicely sensitive to the covariant 
structure of the isovector fields. At each beam energy we see an increase of the ir~ / 7r + and K / K + yield 
ratios with the models NL — > DDF — > NLp — > NLpS. The effect is larger for the K°/K + compared 
to the 7r~/7r + ratio. This is due to the subthreshold production and to the fact that the isospin effect 
enters twice in the two-step production of kaons, see [16] . Interestingly the Iso-EoS effect for pions is 
increasing at lower energies, when approaching the production threshold. 

We have to note that in a previous study of kaon production in excited nuclear matter the dependence 
of the K° / K + yield ratio on the effective isovector interaction appears much larger (see Fig.8 of ref.|15j). 
The point is that in the non-equilibrium case of a heavy ion collision the asymmetry of the source where 
kaons are produced is in fact reduced by the n — > p "transformation" , due to the favored nn pA~ 
processes. This effect is almost absent at equilibrium due to the inverse transitions, see Fig. 3 of ref.|15j. 
Moreover in infinite nuclear matter even the fast neutron emission is not present. This result clearly 
shows that chemical equilibrium models can lead to uncorrect results when used for transient states of 
an open system. 

7 On the Transition to a Mixed Hadron- Quark Phase at High 
Baryon and Isospin Density 

The possibility of the transition to a mixed hadron-quark phase, at high baryon and isospin density, is 
finally suggested. Some signatures could come from an expected "neutron trapping" effect. 

In order to check the possibility of observing some precursor signals of a new physics even in collisions 
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Figure 9: 238 U + 238 U , 1 AGeV, semicentral. Correlation between density, temperature (black values), 
momentum thermalization (3-D plots), inside a cubic cell, 2.5 fm wide, in the center of mass of the 
system. 



of stable nuclei at intermediate energies we have performed some event simulations for the collision of 
very heavy, neutron-rich, elements. We have chosen the reaction 238 U + 238 U (average proton fraction 
Z/A = 0.39) at 1 AGeV and semicentral impact parameter b = 7 fm just to increase the neutron 
excess in the interacting region. In Fig. [9] we report the evolution of momentum distribution and 
baryon density in a space cell located in the cm. of the system. We see that after about 10 fm/c 
a local equilibration is achieved. We have a unique Fermi distribution and from a simple fit we can 
evaluate the local temperature (black numbers in MeV). We note that a rather exotic nuclear matter is 
formed in a transient time of the order of 10 fm/c, with baryon density around 3 — 4po, temperature 
50 — 60 MeV, energy density 500 MeV fm~ 3 and proton fraction between 0.35 and 0.40, likely inside 
the estimated mixed phase region. 

In fact we can study the isospin dependence of the transition densities [17]. The structure of the 
mixed phase is obtained by imposing the Gibbs conditions [62] for chemical potentials and pressure and 
by requiring the conservation of the total baryon and isospin densities 



where x is the fraction of quark matter in the mixed phase. 

In this way we get the binodal surface which gives the phase coexistence region in the (T,pg,p 3 ) 
space. For a fixed value of the conserved charge p 3 we will study the boundaries of the mixed phase region 
in the (T, p#) plane. In the hadronic phase the charge chemical potential is given by p 3 = 2E sym (pB)-^ ■ 
Thus, we expect critical densities rather sensitive to the isovector channel in the hadronic EoS. 

In Fig. [10] we show the crossing density p cr separating nuclear matter from the quark-nucleon mixed 
phase, as a function of the proton fraction Z/A. We can see the effect of the 5-coupling towards an 
earlier crossing due to the larger symmetry repulsion at high baryon densities. In the same figure we 
report the paths in the (p,Z/A) plane followed in the cm. region during the collision of the n-rich 
132 Sn+ 132 Sn system, at different energies. At 300 AMeV we are just reaching the border of the mixed 
phase, and we are well inside it at 1 AGeV. 




P iH \T,H { B H h = P {Q \T,pW), 



Pb = (1 -X)Pb + XPb, 
p 3 = (1 - X )pf + XPs , 



(11) 
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Z/A 

Figure 10: Variation of the transition density with proton fraction for various hadronic EoS parame- 
terizations. Dotted line: GM3 RM F-model [6T] : dashed line: NLp ; solid line: NLpd~ . For the quark 
EoS: MIT bag model with 5 1 / 4 =150 MeV. The points represent the path followed in the interaction 
zone during a semi-central 132 Sn+ 132 Sn collision at 1 AGeV (circles) and at 300 AMeV (crosses). 

« H =(p n -p p )/p B a Q =3(p d -p u )/(p d +p u )=(p d - Pii )/p B 
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Figure 11: EoS of Simmetric/Neutron Matter: Hadron (NLp), solid lines, vs. Quark (MIT-Bag), 
dashed lines. cxh,q represent the isospin asymmetry parameters respectively of the hadron, quark matter: 
a H,Q = 0, Symmetric Matter; cuh,q — 1; Neutron Matter. 




We can expect a neutron trapping effect, supported by statistical fluctuations as well as by a sym- 
metry energy difference in the two phases. In fact while in the hadron phase we have a large neutron 
potential repulsion (in particular in the NLpS case), in the quark phase we only have the much smaller 
kinetic contribution. Observables related to such neutron "trapping" could be an inversion in the 
trend of the formation of neutron rich fragments and/or of the 7r~ /tc + , K°/K + yield ratios for reaction 
products coming from high density regions, i.e. with large transverse momenta. 

Isospin in Effective Partonic Models 

From the above discussion it appears extremely important to include the Isospin degree of freedom in 
any effective approach to the QCD dynamics. This can be easily performed in a two-flavor NJL model 
[63] where the isospin asymmetry can be included in a flavor-mixing picture [M] via a Gap Equation 
like Mi = mi — 4Gi$j — 4G 2 &j, i ^ j, (u, d) , where the =< uu >,< dd > are the two (negative) 
condensates and m u ^ = m the (equal) current masses. Introducing explicitily a flavor mixing, i.e. the 
dependence of the constituent mass of a given flavor to both condensate, via G\ — (1 — a)Go, G 2 = aG® 
we have the coupled equations 

M u = m- 4G $ M - 4aG ($« - &d), 
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Figure 12: NJL results with flavor mixing a = 0.2 and isospin asymmetry (N — Z)/A = 0.4: costituent 
u- d- masses vs.baryon density. 



For a = 1/2 we have back the usual NJL (M u = Md), while small/large mixing is for a =^> 0/a =^> 1 
respectively. 

In neutron rich matter | $^ | decreases more rapidly due to the larger pd and so ($ u — $d) < 0. 
In the "realistic" small mixing case, see also [64, 65J, we will get a definite M u > Md splitting at high 
baryon density (before the chiral restoration). This expectation is nicely confirmed by a full calculation 
[66] of the coupled gap equations with standard parameters (same as in ref.JM]). The results are shown 
in FigfU 

All that represents a more fundamental confirmation of the m* > m* choice in the hadron phase, as 
suggested by the effective QHD model with the isovector scalar S coupling, see before and [12]. However 
this can represent just a very first step towards a more complete treatment of isovector interactions in 
effective partonic models, of large interest for the discussion of the phase transition at high densities. 
We can easily see that the mass splitting effect is not changing much the symmetry energy in the quark 
phase addressed before. However confinement is still missing in these mean field models. Stimulating 
new perspectives are open. 

8 Perspectives 

We have shown that violent collisions of n-rich heavy ions from low to relativistic energies can bring new 
information on the isovector part of the in-medium interaction, qualitatively different from equilibrium 
EoS properties. We have presented quantitative results in a wide range of beam energies. At low 
energies we see isospin effects on the dissipation in fusion and deep inelastic collisions, at Fermi and 
Intermediate energies the Iso-EoS sensitivity of the isospin dynamics in fragment reactions and in 
collective flows. 

We have shown that meson production in n-rich heavy ions collisions at intermediate energies can 
bring new information on the isovector part of the in-medium interaction at high baryon densities. Im- 
portant non-equilibrium effects for particle production are stressed. Finally the possibility of observing 
precursor signals of the phase transition to a mixed hadron-quark matter at high baryon density is 
suggested. 

10 years after the first studies of Iso-EoS effects on reactions [3] the symmetry energy case is still 
open. The results presented in this "Isospin Journey" appear very promising for the possibility of 
exciting new data from dissipative collisions with radioactive beams. 
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Abstract 

Correlations in the nuclear wave-function beyond the mean-field or Hartree-Fock approxima- 
tion are very important to describe basic properties of nuclear structure. Various approaches to 
account for such correlations are described and compared to each other. This includes the hole- 
line expansion, the coupled cluster or "exponential S" approach, the self-consistent evaluation of 
Greens functions, variational approaches using correlated basis functions and recent developments 
employing quantum Monte-Carlo techniques. Details of these correlations are explored and their 
sensitivity to the underlying nucleon-nucleon interaction. Special attention is paid to the attempts 
to investigate these correlations in exclusive nucleon knock-out experiments induced by electron 
scattering. Another important issue of nuclear structure physics is the role of relativistic effects 
as contained in phenomenological mean field models. The sensitivity of various nuclear structure 
observables on these relativistic features are investigated. The report includes the discussion of 
nuclear matter as well as finite nuclei. 



1 Introduction 

One of the central challenges of theoretical nuclear physics is the attempt to describe the basic prop- 
erties of nuclear systems in terms of a realistic nucleon-nucleon (NN) interaction. Such an attempt 
typically contains two major steps. In the first step one has to consider a specific model for the NN 
interaction. This could be a model which is inspired by the quantum-chromo-dynamics[?], a meson- 
exchange or One-Boson- Exchange model [?, ?] or a purely phenomenological ansatz in terms of two-body 
spin-isospin operators multiplied by local potential functions [?, ?]. Such models are considered as a 
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